A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

08(08.05.2; 08.11.1; 08.09.2 HIP 60350; 10.15.2 NGC 3603) 



ASTRONOMY 

AND 
ASTROPHYSICS 



Origin and possible birthplace of the extreme runaway star HIP 
60350 

P. Tenjes 1 2 , J. Einasto 2 , H. M. Maitzen 3 , and H. Zinnecker 4 

1 Institute of Theoretical Physics, Tartu University, Tahe 4, Tartu 51010, Estonia 

2 Tartu Observatory, 61602 Toravere, Estonia 

3 Institut fur Astronomie der Universitat Wien, Turkenschanzstrasse 17, A-l 180 Wien, Austria 

4 Astrophysikalisches Institut Potsdam, An der Sternwarte 16, D- 14482 Potsdam, Germany 

Received 28 September 2000 / accepted date 



Abstract. Using the recently determined spatial velocity com- 
ponents of the extreme runaway star HIP 60350 and a gravi- 
tation potential model of the Galaxy, we integrate the orbit of 
HIP 60350 back to the plane of the Galaxy. In this way, a pos- 
sible location of the formation of the star is determined. We 
estimate the uncertainty of the result due to the uncertainties 
of the gravitational potential model and the errors in the spatial 
velocity components. The place of birth lies (within the errors) 
near the position of the open cluster NGC 3603. However, the 
ejection event which occured about 20 Myr ago is in contradic- 
tion with the cluster mean age of 3 — 4 Myr. We suggest that it 
occured at an earlier phase in sequential star formation in that 
region. We discuss also ejection mechanisms. Due to the rather 
high mass of the star (about 5 M ), the most probable model is 
that of dynamical ejection. 

Key words: stars: formation - stars: kinematics - stars: indi- 
vidual HIP 60350 - open clusters and associations: NGC 3603 



1. Introduction 

Recent Hipparcos proper motion measurements for the star HIP 
60350, together with its radial velocity, allowed the determina- 



tion of its spatial velocity components U 



V 



183 kms" 1 and W 



-130 kms~ 



+352 kms , 
giving for the 



1998) 



total velocity the value v = 417 km s~ (Maitzen et al 
These results suggest the star is a rare extreme runaway star. 
The origin of such a high velocity is not clear, especially given 
the rather high mass of the star, M ~ 5 M . Thus, we have 
chosen to study the possible origin of the velocity: first the 
location of birth, and then the mechanisms allowing the ac- 
quisition of such a high velocity. Maitzen et al. (1998) calcu- 
lated the distance of HIP 60350 from the Sun, about 3.5 kpc. 
Together with its equatorial coordinates, this gives us galacto- 
centric cylindrical coordinates of HIP 60350, R — 9.25 kpc, 
6 = 3.3°, z = 3.4 kpc (zero point of the azimuthal angle is 
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in the direction to the Sun). Hence, both the velocities and the 
coordinates of the star are known. 

The star HIP 60350 is of spectral type B4-5V. The high ve- 
locity carries it far from the galactic plane, towards the outer 
halo. For this reason, the study of its orbit is interesting in re- 
lation to another problem, namely the origin of young stars ob- 
served in galactic halos. Recent surveys of the Milky Way high 
latitude (halo) reg ions i dentified several doze n young B-type 
stars (Saffer et al. 



1997 



Lindblad et al. 



1997 



Rolleston et al. 



1999b - Within 0.7 kpc from the Sun, Hoogerwerf et al. ( |200C| ) 



identified 56 runaway stars. Slightly more uncertain data indi- 
cate t hat th ere are also yo ung st ars in the halo of M 3 1 (Hambly 



1995 



Smoker et al. 2000) 



et al. 

The origin of these stars remains unclear, as well as their 
place of birth. Two main theories exist. According to the first 
theory these stars are real halo stars, i.e. they were born at high 
galactic latitudes as a result of col lision s between high-velocity 
clouds (e.g. Dyson & Hartquist 1983 ). However, for several 
stars of this kind their radial velocities indicate that they are 
moving away from the Galactic plane with high speed. Thus 
the second theory is based on the hypothesis that these stars 
are ejected with high velocities from the disk by some mech- 
anism. Two main mechanisms have been considered so far - 



ejection via supernova explosion (Blaauw 1961 ), and dynami- 
cal ejection via close stellar encounters in star clusters (see e.g. 
Leonard & Duncan 1988). Unfortunately, spatial velocities of 



these stars are poorly known and even if known they may in- 
clude systematic errors (Rolleston et al. 1999j ). 

For this reason, the study of HIP 60350 provides a good op- 
portunity to determine the possible place of birth. In an earlier 
study, Maitzen et al. (1998), using simple dynamical consider- 



ations, estimated the time needed for the star to move from the 
galactic plane to its present position (about 20 Myrs), and spec- 
ulated that the possible place of birth could be the spiral arm 
-II in the 4th quadrant. In the present study we use a detailed 
gravitational potential model of the Milky Way, integrate nu- 
merically the velocity components and coordinates backwards, 
and determine the place of birth in the galactic plane. We es- 
timate also the uncertainties of the result and look for nearby 
star clusters. 
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2. Gravitational potential of the Galaxy 

Numerical calculation of a stellar orbit presumes the knowl- 
edge of the gravitational potential of the galaxy. In the present 
paper, our main interest is not the long-term parameters of an 
orbit, i.e. the behaviour of the orbit in the far future. This kind 
of analysis depends essentially on the mass distribution of the 
Galaxy at large galactocentric distances, particularly on the dis- 
tribution model of dark matter. Rather, our main aim here is to 
find the possible place of birth of the star. For this, it is impor- 
tant to know the mass distribution of the Milky Way mainly 
at intermediate distances. For that region, an accurate rotation 
curve of the Galaxy exists (construction of a mass distribu- 
tion model of the Galaxy includes, beside the rotation curve, 
the knowledge of additional parameters, e.g. the local Oort 
constant s, the solar distance and the circular velocity etc., see 
Einast o (1979 1 or a more recent review by Dehnen & Binney 
( 1998 )). The circular velocity is directly related to the radial 
derivative of the gravitational potential, which is contained in 
the equations of motion of a star. Several models exist which 
describe the gravitational potential of the Galaxy. In order to 
calculate the orbit of the star HIP 60350, we need to know the 
radial and vertical derivatives of the potential also outside the 
galactic plane. Hence, it is better to use models which take into 
account the real ellipticity of the Milky Way mass distribution, 
including the individual flatness of different stellar populations. 
For this reason we decided to use the model of the Galaxy pro- 
posed by Haud & Einasto ( 1989 ). A recent model by Dehnen & 
Binney (1998) also takes into account the real flatness of stel- 



lar populations but their algorithm for gravitational potential 
calculation is less convenient for our purpose. 

In our model, the Galaxy is given as a superposition 
of different subsystems. Each subsystem represents a cer- 
tain stellar/gas population with corresponding density distri- 
bution, chemical composition and kinematical characteristics. 
The density distribution of each component is approximated by 
an inhomogeneous ellipsoid of rotational symmetry with con- 
stant axial ratio e (Einasto 1972, Einasto & Haud 1989] ). The 
spatial density of visible populations is described by the law 

p(a)=p(0)exp[-(a/a c )V N ], (1) 

where p(0) — hMj (47reag) is the central density, a = 
\/i? 2 + z 2 /e 2 is the distance along the major axis, a c = kao 
is the core radius (ao is the harmonic mean radius), h and k 
are normalizing parameters, depending on the parameter N, 
which allows to vary the density behaviour with a. The def- 
inition of normalizing parameters and their calculation is de- 



scribed in Tenjes et al. (1994 ), appendix B. For the disk and the 



flat components we use the density distribution in the following 
form 



P{°) = p+{a) - p-(a), 



(2) 



where subindices "+" and "— " denote density distributions (1) 
of components with positive and negative masses respectively. 
In this way, we obtain density distributions with a central den- 
sity depression. If we demand that the density be zero at a = 



and positive elsewhere, the following relations must hold be- 
tween the parameters of components p + and p_: ao- = «ao+, 
M- — — k 2 M + , e_ = e+/n, where k < 1 is a parame- 
ter which determines the relative size of the hole in the cen- 
tre of the disk. To avoid negative densities of the population, 
the structural parameters of components N+ and iV_ must be 



equal (see Einasto et al. 1980) 



The dark matter distribution is represented by a modified 
isothermal law 



p{a) = J /*0)([1 + (£) 

lo 



_a_ N i2l-l/2 



[l + (g)2]-l/2)2 a < fl 

a > a . 

(3) 

Here a is the outer cutoff radius of the isothermal sphere. 

The modeling procedure and model parameters are given 
in the original paper and it is not necessary to repeat them 
here. Knowing the density distribution formula of the compo- 
nents and the parameters of the populations we can calculate 
the gravitational potential at every point (R, z). 

3. Numerical orbit calculation 

The equations of mo tion in cylindrical coordinates are (see e.g. 
Binney & Tremaine |1987[ Sect. 3.1) 



(4) 



L z = Rv^p = const, 

8R + R 3 ' 

oz 

where <I> is the gravitational potential not depending on the az- 
imuthal coordinate ip. The gravitational potential derivatives 
for inhomogeneous ellipsoidal mass distribution with constant 
ellipticity are 

d$(R,z) GhM /- aresin ( e ) 

— R- — / p (a) sin x dx 



dR 

d$(R,z) 
dz 

where a 2 = — 



(ea ) 3 

GhM 
(ea ) 3 



(5a) 



arcsin(e) 



p* (a) tan 2 x dx (6a) 



\ COS z X j v 



e 2 is the eccen- 



tricity and p*(a) = exp [a/lkao)] 1 ^^ . Here p*(a) is our 
density distribution (1). 

For a spherical dark matter distribution the gravitational po- 
tential derivatives are 

d$(R,z) GhM f 1 , 2 

= R — — I p*{a)u z du, (56) 



dR 

d$(R, z) 
dz 



GhM 



p*{a) u 2 du, 



(66) 



*o Jo 

where a 2 = u 2 (R 2 + z 2 ) and p*(a) is again the density distri- 
bution (3) without the factor p(0). 

We calculated integrals (5) and (6) numerically using the 
Gaussian quadrature formula with 40 points. The equations of 
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Fig. 1. Test orbit with (v R , v^,v z ) = (0, 217, 5). 

motion (4) were solved numerically using the 4th order Runge- 
Kutta method. 

The correctness of the calculated orbits was tested in sev- 
eral ways. The most simple test is to calculate the circular ve- 
locity at some distance (e.g. at solar distance) with the known 
formula V 2 — Rd<& /dR, and then integrate the orbit numer- 
ically with these initial data. Doing so we found that the orbit 
remains circular with a precision of 0. 1 percent. 

Next we calculated vertical oscillations relative to the 
galactic plane at solar distances with v z = 5 km s _1 at z = 0. 
The orbit is presented in Fig. [l]. We see that the period of 
vertical oscillations remains constant, as it must be, and is 
T = 0.0896 Gyrs. For small oscillations this period is simply 
related to the local mass density: 

7T A 2 - B 2 



GT 2 



2nG 



Here A = 14.8 km s 1 kpc~ 1 and B = -12.4 km s _1 kpc" 



are Oort's constants (Feast & Whitelock 1997 1, giving the value 
p = 0.0890 M©pc~ 3 . The direct calculation from the initial 
parameters of the galactic model yields as the local mass den- 
sity 0.0884 M Q pc -3 , in very good agreement with the previ- 
ous value. 

Finally, we calculated a typical orbit with {vr^v^^Vz) — 
(50, 250, 50). This orbit in the comoving meridional plane is 
presented in Fig. ^ . We see that the orbit remains well confined 
within a region limited by the integrals of motion (see Kuzmin 



1956,OHongren 1962) 



After these tests we began to calculate the orbit of HIP 
60350. 



4. Orbit of the star HIP 60350 

Knowing the present phase space coordinates for HIP 60350 
we calculated its orbit back, and found that the star was in the 
galactic plane 20.4 Myrs ago at R — 8.05 kpc and 6 = -80°. 



Fig. 2. A meridional section of a test orbit with (vr, ity, v z 

(50,250,50). 



The spatial velocity components at that moment were v r — 
-278 km s" 1 , v s = 463 km s"\ v z = 191 km s^ 1 . Subtract- 
ing the galactic rotation component at 8 kpc v# = 220 km s _1 
we obtain that the ejection velocity was v = 416 km s" 1 . The 
projection of the orbit into the plane of the Galaxy is shown in 
Figure |3| by a bold solid line. 

In order to examine possible ejection mechanisms it is 
needed to estimate the uncertainties in our result. First, we es- 
timate errors due to the galactic mass model. An acceptable 
model must be in accordance with measured galactic rotation 
velocities. Thus, the mass of the galactic disk (which is the 
most essential parameter in our case) must give rotation ve- 
locities within the measured velocity errors. The rotation curve 
at distances R ~ 8 — 10 kpc from the Galactic centre is known 



with errors ±10 km s (Fish & Tremaine 1990, Binney & 



Merrifield 1998). These errors allow us to vary the disk mass 
within the limits ±0.9- 10 10 M . Orbit calculations for galactic 
models with disk masses 5.9 • 10 10 M and 7.8 • 10 10 M© show 
that the position of galactic plane crossing has an uncertainty 
of only ±0.1 kpc, which is rather small. 

Second, we studied uncertainties due to the errors in the ob- 
served velocity components of the star, quoted by Maitzen et al. 
(1998) to be 15 percent. To the ^-component velocity we must 
add the circular velocity of LSR (220 ± 10 km s _1 ). Thus the 
uncertainties of the velocity components are (±50, ±30, ±20). 
To estimate the influence of velocity uncertainties of each com- 
ponent on the position of the star at the galactic plane crossing, 
we calculated orbits of the star for a number of observed veloc- 
ities including estimated errors. We find that the resulting over- 
all uncertainty of the plane crossing lies approximately within 
an ellipse with semiaxes 1.1 kpc and 0.7 kpc. The long axis 
of the ellipse is oriented approximately in the direction of the 
stellar orbit (Fig. ||, dashed ellipse). Hence we must search for 
the possible birthplace of HIP 60350 within this region. 
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Fig. 3. Motion of clouds Georgelin 56, 57, 60 and 62 (=NGC 
3603), Grabelsky 19, 22, and 27, the Sun and the star HIP 
60350 at the Galactic plane during 20.4 Myrs. For HIP 60350 
the orbit projected to the Galactic plane is given. 



Because our star is young, we are interested in the pres- 
ence of young star clusters and associated interstellar gas re- 
gions. According to a catalog of HII regions of the Milky Way 



by Georgelin & Georgelin (1976 1, potential candidates could 



be objects Nos. 56, 57, 60 and 62. In addition, we looked for 
promising candid ates in the molecular gas survey by Grabel- 
sky et al. ( 1988 >; clouds Nos. 19, 22 and 27 are candidates. 
From these lists only clouds No. 60 and 62 are known to con- 
tain stars at present. Moreover, cloud No. 62 corresponds to 
the well-known star cluster NGC 3603; thus, the distance to 
this cluster is known with high accuracy. The distances to other 
clouds without known associated stars are estimated by the au- 
thors of corresponding papers using indirect data and are less 
certain (distances to the cloud s Nos 56 and 57 are estimated by 
Georgelin & Georgelin ( 1976) on the basis of galactic northern 
and southern rotation models, distances to the clouds Nos 19, 



22 and 27 are estimated by Grabelsky et al. (1988) on the basis 
of radius-line-width relation; in both papers the authors used 
the old values of galactic constants (10 kpc and 250 km s _1 ), 
thus we transformed the clouds distances to the new values 
of galactic constants recommended by the IAU, see Kerr & 



Lynden-Bell (1986)). As the star HIP 60350 started from the 



galactic plane about 20 Myr ago we must transfer the positions 
of the selected clouds back by that time assuming circular or- 
bits at their galactocentric distances. The results are presented 
in Fig. [3} 

The positions of young objects in the galactic plane may 
be influenced by their peculiar velocity component in the R- 
direction. Peculiar velocities of young star clusters are typically 



< 15 kms" (see e.g. Nezhinskij et al. 1995), giving during 
a 20 Myr orbit a maximum correction of the position up to 
±0.3 kpc. Thus we see that clouds Nos. 60 and 62, which are 
known to contain stars, as well as molecular clouds, may lie 
inside the boundary of the permitted region of the birthplace of 
our star. 

Recent observations (Brandl et al. 



1999, De Pree et al. 



1999) allowed us to study in detail the star cluster NGC 3603. 



According to these observations this cluster lies at a distance 
6.1 ± 0.6 kpc from the Sun. NGC 3603 contains a lot of young 
and massive stars and is still in the stage of star formation. Its 
position corresponds to the birthplace of HIP 60350 nearly per- 
fectly. Related to the cluster HII region is one of the largest HII 
regions in the Galaxy. The initial mass function (IMF) of the 



cluster extends up to 120 M Q (Drissen et al. 1995 1. 



5. Discussion 

In the present work we integrated the orbit of the extreme run- 
away star HIP 60350 back to the galactic plane. Taking into 
account possible errors in the velocity measurements and the 
galactic model we have found the probable birthplace of this 
star. 

OB runaway stars can be produced either by supernova ex- 
plosions in massive close binaries or by close dynamical en- 
counters in dense star clusters. First we discuss the ejection as 
a result of SN explosion. This mechanism was studied first by 
Blaauw (1961) as a symmetric SN explosion where the mo- 
mentum of ejected matter was balanced by the momentum of 
the binary centre-of-mass. However, before the SN explosion 
there must occur some mass transfer between the system com- 
ponents, and thus the amount of ejected matter is less than half 
of the system's total mass. Disruption of a binary system is un- 
likely and the maximum velocity of the centre-of-mass cannot 
be too large. According to Tauris & Bailes (1996), the limit 



for the recoil velocity is 270 km s (see also Nelemans et al. 



1999) 



According to recent studies, supernova explosions can be 
asymmetric. To explain the observed spatial velocity distribu- 
tion of single pulsars, we need to assume asymmetric SN ex- 
plosions, where the newborn neutron star, due to an additional 
kick, attains a mean velocity of 450 km s" 1 in an arbitrary 



direction (Lyne & Lorimer 1994, Hartmann 1997). This ex- 



plosion probably disrupts the system, and thus runaway stars 
should not be presently binary stars. Although an asymmet- 
ric explosion may give a very high velocity to the pulsar (up 
to 1500 km s _1 ), its impact to the companion star is signifi- 
cantly smaller. According to calculations by Tauris & Takens 
( 1998 ) the spatial velocity of the companion star depends on 
several parameters (kick speed, companion star mass etc); usu- 
ally it does not exceed 300 km s~ . Applying these results to 
our 5 Mq star, the runaway velocity is less than 200 km s , 
even for the largest kick speeds. 

Because it is difficult to obtain the ejection velocity of 
420 km s _1 with a SNe scenario, we turn our attention to dy- 
namical interactions. Dynamical interactions have been stud- 
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ied by Leonard & Duncan (1988), who showed that typical 
velocities in interactions extend up to 250 kms -1 . Accord- 
ing to models the most efficient in producing runaway stars 



are binary-binary interactions (Mikkola 1983). The efficiency 
of binary-binary interactions dominates over binary-single in- 
teractions, especially in producing high-velocity runaway stars 



(Leonard & Duncan 1988). However, these are the results of N- 



body experiments and were limited by available computational 



time . For this reason Leonard (1991) performed special numer 



ical experiments with binary-binary interactions to determine 
the maximum possible velocity of ejection. His results were 
presented for different ratios v e j/v esc (v e j - ejection velocity 
of a star, v esc - escape velocity at the surface of ejected star). 

HIP 60350 was ejected from the galactic plane with a spa- 
tial velocity of 420 km s -1 . A 5 M Q star has a radius of 2.6 Kq 
(Tout et al. 1996) and v esc = 840 km s~~ , thus v e j — 0.50u esc . 



According to Leonard (1991, his Table 3) the ejection veloc- 
ity v e j = 0.50w esc cannot be obtained if HIP 60350 inter- 
acted with equal or lower mass stars. Such high ejection ve- 
locities can be obtained by a star in interactions with three 
stars having masses at least 4 times larger, but probably even 
8 times larger. In this case sufficient ejection velocities for the 
least massive companion may appear within about 400 cluster 
crossing times. Hence to produce a high-velocity star like HIP 
60350, the participation of 3 stars with masses ~ 40 M is 
needed. Typical diame ters of young open clusters are 1 — 9 pc 
(Phelps & Janes 1993 ), the line-of-si ght ve locity dispersion is 
about 3 — 4 km s - (see e.g. Rroupa 200C ), thus the resulting 
crossing time is 0.3 — 2 Myr. On the other hand, according to 
photometric measurements, HIP 60350 is clear ly rath er close to 
the ZAMS, and cannot be old (Maitzen et al. |1998| ). It is very 
plausible that the ejecting event took place in an early phase 
of star formation in a very compact cluster, where IMF was 
skewed towards the high-mass end compared with the field star 
IMF (Clarke & Pringle 1992). If this is the case, stars of larger 
mass could indeed participate in the ejection event (MS life- 
times of 20 M Q and 40 M Q stars are 5 Myr and 2 Myr, respec- 
tively). 

Dynamical interactions are more effici ent in massive and 
compact clusters (Leonard & Duncan 1990 ). Interactions (ejec- 
tions, mergers) are also more frequent in clusters with ac- 
tive s tar formation at early e voluti onary stages (Bonnell et al. 



1998 , Portegies Zwart et al. 1999 ). For this reason the cluster 
G 62 = NGC 3603 seems to be a suitable candidate. The IMF 
of the cluster extends up to large masses, making high ejec- 
tion velocities reasonable. When varying the observed veloc- 
ity of HIP 60350 (but remaining within observational errors) 
we calculated that with the present velocity components v r = 
318 km s~ , vg = 423 km s -1 and v z — 138 km s -1 the place 
of birth of the star coincides with the position of G62 (Fig. ^j, 
circle and continuous line, see also Table 1, Model 1). 

However, there exists a serious argument against NGC 
3603 as the birthplace of HIP 60350. The age of NG C 360 3 
is estimated to be only 3 — 4 Myrs (De Pree et al. 1999), 
which is in contradiction to an ejection 20 Myrs ago. Further, 
some time is also needed for interactions - statistically up to 
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Fig. 4. Positions of clouds at the times of star ejection. Corre- 
sponding orbits of the ejected star are labeled as different mod- 
els with parameters given in Table 1 . 

Table 1. Parameters of stellar orbits for individual clouds 



Model Present velocity (km/s) Ejection 



Ejection 







ve 


v z 


time (Myr) 


velocity 


1 


318 


423 


138 


19.9 


405 


2 


307 


428 


130 


20.7 


402 


3 


332 


423 


130 


20.7 


416 


4 


343 


418 


130 


20.5 


420 


5 


397 


403 


117 


21.6 


456 



400t C ross > 100 Myr. This discrepancy may not be critical, be- 
cause star formation in young clusters is a complicated process 



with different stages (see e.g. Elmegreen 2000 and for NGC 



3603, Brandneret al. 1997 1. Thus, we would like not to exclude 
this cluster together with its surrounding region as a candidate 
for the origin of HIP 60350. The age structure of t he clu ster 
deserves special and careful study: Eisenhauer et al. 1998 de- 
rived that the distribution of stellar ages in NGC 3603 is non- 
Gaussian and extends up to 100 Myr. 

Another candidate is the HII region G60, which is also sit- 
uated quite near to the probable ejection place. With present 
velocity components v r = 307 kms -1 , vg — 428 kms -1 
and v z = 130 km s -1 the ejection place will lie close to G60 
(Fig. [|, Table 1, model 2). The HII region G60 is less studied. 
We do not know the detailed structure and age of this region. 
These properties remain suitable for further study. 

According to the study of molecular clouds by Grabelsky 



et al. (1988), near the birthplace of HIP 60350 there are clouds 
of molecular gas No 19, 22 and 27 according to his desig- 
nations (in addition cloud No 17 corresponds to NGC 3603). 



6 



P. Tenjes et al.: Origin and possible birthplace of the extreme runaway star HIP 60350 



400 
300 
200 
100 


-100 
-200 
-300 
-400 




tion in the adjecent parts of the original cloud. These events 
might be similar to the scenario outlined by Preibisch & Zin- 
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Fig. 5. The orbit of the star HIP 60350 for the next 13.5 Gyrs, 
projected to the galactic plane. Circle with centre coordinates 
(0, 0) is situated at the galactic centre and has a radius 10 kpc. 
The x and y axes are directed toward the direction of rotation 
(at Sun position) and away from the centre toward the Sun, 
respectively. 

For these clouds the corresponding models are 3-5, see Ta- 
ble 1. Unfortunately the distances to these clouds are kine- 
matic (they have based on empirical radius-line-width relations 
by Dame et al. ( 198fj| ), see Grabelsky et al.) and are thus un- 



certain. Moreover, these clumps may be bl ended together at 
lower intensities, and SN explosions (Opik 1954 ) and stellar 
winds from massive stars (Carpenter et al. 1995) may disrupt 
the original cloud structure. Molecular clouds have lifetimes of 
10 - 100 Myr (Blitz & William s [2000| , Williams et al. |2000| ) or 
less (Ballesteros-Paredes et al. 1999, Elmegreen 2000), which 
is probably too short to produce high ejection velocities. Thus, 
it is not likely that the referred molecular clouds can be directly 
assigned as birthplaces of HIP 60350. 

At present the place of birth of our star is at I = 298, b = 
and R = 6.9 kpc from the Sun. This region (see also the ellipse 
of errors in Fig. ||) deserves further careful study. According to 
our analysis it is probable that in this region there is an open 
cluster with an age of more than 20 Myr. Finding this kind of 
star cluster is not simple, due to the relatively large distance and 
possible ISM obscuring (unless we are lucky enough to find it 
in a cloud hole). In order to decrease the probable search area 
from where the ejection took place it seems most promising to 
determine more precise velocity measurements for the star. 

Individually identified molecular clouds referred to above 
may simply be remnants of an earlier star formation process. 
Together with NGC 3603 they may form a bigger stellar- 



gaseous complex. Observations by Grabelsky et al. (1988 ) hint 



that near NGC 3603 a more massive and extremely disrupted 
cloud complex exists which refers to past intensive star forma- 
tion. Massive stars from this first star formation stage ejected 
HIP 60350, exploded as supernovae and triggered star forma- 



necker (1999) for Scorpius-Centaurus OB association. 



Observations of different individual OB runaways support 



both the supernova explosion scenario (e.g. Kaper et al. 1997) 
and the cluster ejection model (e.g. Ryans et al. 1999, Mof- 
fat et al. 1998); both scenarios produce runaways (see recent 



detailed analysis about nearby runaways by Hoogerwerf et al. 



2000). The very high velocity runaway HIP 60350 seems to 



been ejected via dynamical ejection from G62 or from G60. 

Finally, we have calculated the future orbit of the star 
HIP 60350, starting from its present position. The total galac- 
tocentric velocity of the star is less than the escape velocity at 
R = 9 kpc, 600 km s , thus it must remain within the po- 
tential well of the Galaxy. The projection of the orbit to the 
galactic plane is shown in Fig. [|. We see that its apogalactic 
distance reaches in the first and second revolution 425 and 554 
kpc at 3.3 and 1 1.4 Gyrs, respectively. Evidently, at these large 
distances, perturbations by other members of the Local Group 
galaxies are important. Most likely the star will remain in the 
common potential well of the whole Local Group. 
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